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Abstract. 

We study the faint end of the HI mass function (HIMF) in order to test the predictions of the 
CDM theory on the number density of objects with small (dark) masses. The neutral hydrogen 
is much better tracer of the underlying mass distribution compared to the luminous matter and 
can be used to test the existence of a population of small galaxies in which the star formation has 
been partially or completely suppressed during cosmic evolution. Due to technical limitations, 
the existing HI surveys are not very sensitive to HI masses below 10* Mq . We designed a blind 
HI survey to be sensitive to objects with small HI masses. The surveyed area is in the Canis 
Venatici groups of galaxies and covers in total ~ 86 deg^ of sky, with observed velocities in the 
range —350 < cz < 1400 km s~^. We detected 69 objects, 22 of them for the first time in HI. All 
new HI detections fall in the lower part of the mass-histogram, confirming our ability to detect 
galaxies with small HI masses. The calculated HIMF is fiat in the faint end regime (a ~ -1), 
different from the steep rise predicted by CDM models. Possible effects of the environment on 
the estimated HIMF parameters are discussed. 
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1. Motivation 

Some works that address the discrepancy between the number of the observed dwarf 
galaxies in the optical and the number of those galaxies predicted by CDM models, 
suggest that small galaxies contain relatively less luminous matter than larger galaxies 
(e.g. Verde, Oh & Jimenez 2002). If that is the fraction of galaxies with dark 

matter haloes of mass below 10^ - 10^° Mq may never form stars and instead they will 
remain dark. Neutral hydrogen, if still within these systems, would be a much better 
tracer of dark galaxies. If one assumes that HI makes up a few percent of the total mass 
of a galaxy, dark galaxies would contain HI in the range 10"^ - 10^ Mq or less. The 
surveys carried out in HI have poor statistics for objects with masses below 10^ M©. 
Moreover, all HI objects detected on extragalactic scales have an optical counterpart and 
have masses above 10^ Mq (Briggs 2004). On the other hand, there is a population of 
recently discovered compact HI clouds (CHVCs) around the Milky Way and M31 (de 
Heij et al. 2002). The CHVCs have HI masses below 10^ Mq and none of them has an 
optical counterpart. The gap in the observed HI masses between 10^ - 10^ M© is due 
to technical limitation of existing surveys. The existing HI information is not sufficient 
in the range of HI masses predicted in these galaxies with little or no stars. There is 
obviously a need for an HI survey in which the galaxies with HI masses below 10^ Mq 
are a significant fraction of all detections. 
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Figure 1. Number of detections per HI mass bin. Tlie sliaded histogram represents the distri- 
bution of the whole sample of 69 detections, the filled histogram represents the distribution of 
22 detections observed for the first time in HI. 

2. The CVn blind survey 

We designed an HI survey which is extremely sensitive for objects with HI masses 
below 10^ Mq. The main difficulty for the HI observations are the long integration 
times. To be able to detect a statistically meaningful sample of objects with low HI 
masses in a reasonable amount of time, we have chosen to observe in a slightly over- 
dense region instead of surveying the field. The Canis Venatici (CVn) groups of galaxies 
are an excellent region for this purpose. The CVn galaxies extend along the line of sight 
up to a distance determined by cz « 1200 km s~^. This region contains a large number 
of the dwarf galaxies, rich in HI, as detected in the previous surveys (in the optical by 
Binggeli et al. 1990, in the HI by Kraan-Korteweg et al. 1999). 

For the survey, we observed 60 x 12 hr using the Westerbork Synthesis Radio Telescope 
(WSRT) in mosaic mode. The observed area on the sky is « 86 deg^. It is covered 
by observing 60 strips of different declination, separated by 15 arcmin. Each strip of 
constant declination has 24 pointings separated by 15 arcmin in Right Ascension. The 
total integration time per pointing was 80.1 min. Only data with —350 < cz < 1400 
km s~^ were used in further data reduction and analysis process. Given that the size of 
the WSRT primary beam is 34 arcmin, we obtained nearly uniform sensitivity over the 
whole observed area. This makes the CVn survey extremely suitable to detect objects 
with small HI masses. 

The data were calibrated and analysed using scripts written by Oosterloo and Kovac 
based on the the MIRIAD (Sault at al. 1995) routines. The output of this process is ~ 
1400 3-dimensional data-cubes, each cube being created by extracting the central part of 
the combined 9 nearest pointings. The spatial resolution of the data-cubes produced is 
~ 30 X 60 arcsec. The velocity resolution is ~ 33 km s~^. The noise in the cubes is just 
below 0.9 mjy Beam^^. The data-cubes were automatically searched for signal, using 
algorithms implemented by Kovac. The performance of the algorithms has been tested 
by inserting simulated objects into the real data-cubes. The HI mass of the detections 
was estimated using Mhi = 236D^ f Sdv[MQ] , where D is the distance to the object in 
Mpc, S is the flux in mJy integrated over the velocity range with significant detection, 
and the width of the channels is dv given in km s~^. 

We used the minimum detectable integrated flux to estimate the maximum distance 
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up to which the object with a certain HI mass can be detected in our survey. Using 
the real and the simulated detections we calculated that S'min = 0.2 Jy km s^^ is the 
minimum detectable integrated flux to classify a detection in our survey as a real object. 
The minimum HI mass which an object must have to be detected in our survey is then 
given by logMHi,min = 4.67 + 21ogZ?[M0] . Objects with HI masses 10^ h^^ Mq can be 
detected up to a distance 4.6 H^q Mpc and with HI masses 10^ H^q Mq up to 14.6 H^q 
Mpc. 

In total, we detect 69 objects with HI emission. The cross-correlation with objects 
from the existing catalogues has been done using the astronomical data bases NED 
://nedwww. ipac.caltech.edu) and HyperLeda (http://leda.univ-lyonl.fr). We found 
that 4 of our HI detections have not been previously catalogued and 1 object does not 
have an optical counterpart. The last detection is most probably a tidal satellite. It is 
important to point out that 22 of our 69 detections have not been previously detected in 
HI. The number distribution of detections per HI mass bin is shown in Figure 1. 

At this moment the largest uncertainty in the estimated HI masses of the detected 
objects is their distance. We obtained the distances by measuring recessional velocities 
and correcting them for the motion of the Sun in the Galaxy, motion of the Galaxy in 
the Local group and for the infall of the Local Group towards Virgo, using the infall 
velocity of the Local Group according to Theureau et al. (1998). 



3. The HI Mass Function 

The HI mass function (HIMF) is defined analogously to the optical luminosity function, 
so one can use identical methods for its derivation. We used the T,{l/Vmax) method 
introduced by Schmidt (1986). This method consists of summing reciprocal values of the 
maximum surveyed volume in which the objects of a certain HI mass can be detected. 




Figure 2. Left panel: The HIMF from our survey compared to the HIMFs derived from three 
recent blind HI surveys. Our function is re-normalized to match the high mass end of the HIMF 
determined from the plotted surveys. Right panel: Comparison of the HIMF from our survey 
to the HIMF derived for regions of three different local matter densities. The higher values of n 
correspond to regions of higher density (Springob et al. 2005). 

We fitted the binned data with a Schechter function. The best-fit values are a = —1.07, 
log(M*/MQ) = 9.73 and </) = 0.129ft.fo Mpc^^. The calculated HIMF is presented in 
Figure 2. In the left panel of Figure 2, the HIMF from this survey is re-normalized to 
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match the the high mass end of the HIMF derived from three recent bhnd HI surveys. 
The real number density of the HI objects is 20 x higher than plotted. For the surveyed 
volume, the number density of the objects with HI masses spanning ^ 2 decades (10^ 
- 10^ in Mq) is approximately constant. This result is not consistent with the slopes 
of the HIMF as found by Zwaan et al (2005) and Rosenberg & Schneider (2002). The 
disagreement is even larger when comparing this slope to the analytically derived slope 
~ -1.8 of the low mass end of the CDM mass function (Press & Schechter 1974). 

Recently, a lot of effort has been spent to explain the effect of environment on the 
properties of galaxies and their distributions. The HI surveys cover a much smaller area 
than the optical surveys. The only way to study the environmental effects on the HIMF 
is by direct comparison of the estimated HIMF from the surveys carried out in different 
environments (clusters, groups, voids). The main conclusion from this comparison is that 
HI galaxies detected in a denser environment have HIMF characterised with the flatter 
end compared to the HIMF derived from the field. Springob et al. (2005) investigated 
the HIMF of optically selected HI galaxies. Their sample was large enough to investigate 
the dependence of the HIMF on the environment. Since their subsamplcs of galaxies in 
regions of different densities had the same fractions of morphological types, the difference 
in the derived HIMF must be due to the density of the environment. Although the results 
obtained had statistically small differences, the conclusion was that the HIMF in the 
higher- density region shows a flattening in the faint end and lower values of M». Thus, 
the result of the HIMF from our survey is in agreement with this flnding (Figure 2, right 
panel). 
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Discussion 

Karachentsev: About 50 % of the area surveyed by you lies within the Arecibo zone. 
This area will be surveyed soon with much deeper detection limit. It would be reasonable 
to shift your strip of survey at higher declination, i.e. outside the Arecibo zone. 



